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(Dated: March 23, 2007) 

Abstract 

We have studied semileptonic B decay to the exclusive charmless states ir, p/cj, rj and rj' using 
the full 15.5 fb _1 CLEO T(45) sample, with measurements performed in subregions of phase 
space to minimize dependence on a priori knowledge of the form factors involved. We find total 
branching fractions B{B° -► -K~i + v) = (1.37 ± 0.15 stat ± 0.11 sys ) x 10~ 4 and B(B° -► p~t+v) = 
(2.93 ± 0.37 stat ± 0.37 sys ) x 10~ 4 . We find evidence for B+ -► rf£+v, with B(B + -► rf£ + v) = 
(2.66 ± ±0.80 sta t ± 0.56 sys ) x 10~ 4 and 1.20 x 10~ 4 < B{B+ -► rf£+v) < 4.46 x 10~ 4 (90% CL). We 
also limit B(B + — > rj£ + u) < 1.01 x 10 -4 (90% CL). By combining our .B — > 7r£^ information with 
unquenched lattice calculations, we find \V u b\ = (3.6 ± 0.4 stat ± 0.2 sys ^Q^ th ) x 10~ 3 . 
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I. INTRODUCTION 



V U b remains one of the most poorly constrained parameters of the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix [T]. Its magnitude, \V u b\, plays a central role in testing the consis- 
tency of the CKM matrix in B and K meson decay processes. Inconsistency would signal 
existence of new classes of fundamental particles or forces. A precise determination of \V u b\ 
has been the subject of considerable theoretical and experimental effort for well over a 
decade, and remains one of the highest priorities of flavor physics. 

Measurements of exclusive charmless semileptonic decays B — > X u £v [21 El HI El El |8] 
provide a route to the determination of \V u b\ with experimental and theoretical uncertainties 
complementary to the current inclusive techniques (9]. The primary challenge facing all 
measurements of semileptonic b — > u£v decay is separation of the signal from the much larger 
b — > c£v background. Exclusive b — > u£v decays provide unique kinematic constraints for 
this purpose, particularly if a reliable estimate of the neutrino four-momentum is obtained. 

In addition to the experimental challenges facing the exclusive measurements, extraction 
of | V u b | requires knowledge of the form factors that govern the dynamics of these decays. 
The form factors contribute both to the variation of rate, or shape, over phase space and 
to the overall rate via the form factor normalization. The form factors are inherently non- 
perturbative QCD quantities and considerable effort has been devoted to their calculation. 
A major recent theoretical advance has come in the form of predictions from unquenched 
lattice QCD (LQCD) calculations [10J, which provide predictions with uncertainties claimed 
to be at the 10% level for B — > n£v. Light Cone Sum Rules (LCSR) techniques also provide 
predictions for the B — > tc£v [TT] and B — > piv form factors [12J. 

For the extraction of \V U b\, we consider semileptonic decays to the pseudoscalar final 
states ir~£ + v and tt°£ + v, for which we have the unquenched lattice QCD predictions. We 
also investigate the decays to vector final states p~£ + v, p°£ + v, and u£ + v. An accurate 
determination of the decay rate to vector final states is necessary for quantifying cross feed 
into the B — > ir£v rate that is ultimately used to extract \V u b\. 

In addition, we search for B + — > i]£ + v and B + — > r]'£ + v. Because semileptonic decays 
involve only a single hadronic current, they additionally offer a probe of the QCD phenomena 
underlying fully hadronic decay. In particular, the B — > r]£ + v and B — > r]'£ + v decays probe 
the singlet or QCD anomaly component of the rf meson [13]. This component may be 
responsible for the unexpectedly high branching fractions observed in B —>■ 1]'X S decays 
[Ml [15l [16]. A measurement of B(B + — > r]£ + v) was included in previous CLEO studies [I], 
and in 2006, BaBar released preliminary upper limits for both B — ► r]£v and B —>■ rj'£v 
branching fractions |17j . 

This paper describes the CLEO studies of these seven exclusive b — > u£v decay modes. 
The excellent hermeticity of the CLEO detector and the symmetric e + e~ beam configuration 
of the Cornell Electron Storage Ring (CESR) allows for precise determination of the neutrino 
four-momentum. In contrast to asymmetric B factories, the lack of a center of mass boost 
along the beam direction enhances the fraction of T(AS) — > BB events in which all particles 
are contained in the detector fiducial region. This allows more stringent selection criteria to 
be placed on the reconstructed neutrino that improve neutrino resolution without degrading 
efficiency This ultimately results in a cleaner separation between signal and background 
that enables measurements of exclusive semileptonic rates at CLEO to be competitive with 
higher statistics analyses carried out at asymmetric B factories. 

We study all seven decay modes simultaneously, extracting independent measurements 
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in subregions of phase space in order to minimize a priori dependence on form factor shape 
predictions. While the general method is similar to the previous CLEO study [I], the 
current analysis further reduces the sensitivity to the decay form factors by decreasing the 
lower bound on the charged-lepton momentum requirement for the vector modes from 1.5 
GeV/c to 1.0 (1.2) GeV/c for electrons (muons) and, in addition, measuring partial rates 
for vector modes in bins of lepton decay angle. Furthermore, our measurements of the rates 
over the plv phase space allow us to test the validity of form factor shape calculations in 
that mode. 

To increase signal efficiency, we have also modified the phase-space binning in the region 
of high meson recoil momentum, where continuum e + e~ — > qq backgrounds contribute. 

This analysis utilizes a total of 15.4 x 10 6 BB decays obtained at the T(4S), a relative 
increase of 60% from the previous CLEO analysis [4J. The results presented here supersede 
those of the previous analysis. 

The paper is organized as follows. Sections [TT] and |III| discuss the role of form factors for 
extraction of \V u b\ and the relation of the QCD anomaly to the B — > r\tv and B — > rj Iv 



branching fractions. Section IV describes our reconstruction technique. Section |V| describes 
our fitting method, fit components, and fit results. Section IVI] describes the systematic 



uncertainties. Section |VII| interprets the results. Throughout the paper, charge conjugate 
modes are implied. A more detailed description of the analysis technique can be found in 
Reference ITT 
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FIG. 1: Predictions for dT(B — > Trlv)/dq 2 for a variety form factor calculation techniques |10 | lllj , 
120] [21] illustrate the range of variation of the predicted (/ 2 -dependence. 
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FIG. 2: Predictions for dT(B piv)jdq 2 (left) and dT(B p£v)/dcos6 W i (right) for a variety of 
form factor calculation techniques |12[ [2T| [22) [23] illustrate the range of variation of the predicted 
q 2 and cos Oyyi dependence. 



II. EXCLUSIVE CHARMLESS SEMILEPTONIC DECAYS 



For B — > V u £u, where V u is a charmless vector meson, the partial width is 
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Here, k 
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/c is the V u momentum, q is the mass squared of the virtual W (W*), Cq 



(Se) is the cosine (sine) of the angle 6wi between the charged lepton in the W* rest frame 
and the W* in the B rest frame, and C± — 1 ± Cg. H± and H are the magnitudes of the 
W helicity amplitudes, which can be expressed in the massless lepton limit in terms of three 
g 2 -dependent form factors [TH], g, a + , and /, as 



\H ± \ 2 = [f(q 2 ) T 2M B kg(q 2 )}- 
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For a final state pseudoscalar meson P u , H± 
factor f + (q 2 ): 



dT (B 



dq 2 



ub 



0, and the rate depends on a single form 

(4) 
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The structure of the differential decay rates allows us to draw some general conclusions 
regarding the properties of the semileptonic decay modes studied here. For the p{u)iv 
transitions, the left-handed, V — A, nature of the charged current at the quark level man- 
ifests itself at the hadronic level as > \H + \. The H_ contribution is also expected to 
dominate the H contribution, leading to a forward-peaked distribution for cos 9wi- The 
pseudoscalar modes exhibit a sin 2 9wi dependence, independent of the form factor since the 
rate is independent of 9wi- They also depend on an extra factor of p 2 x u , which suppresses 
the rate near g 2 ^ {px u = 0). Taken together, these two effects give the pseudoscalar modes 
a softer charged lepton momentum spectrum than the vector modes. 

Recent results for f+(q 2 ) from lattice calculations with dynamical quarks [10J provide 
a marked theoretical advance in form factor calculations. However, significant theoretical 
uncertainties still remain in the q 2 dependence of the form factors, particularly in the non-7r£^ 
decay modes. Figures [T] and [2] illustrate the variation of form factor predictions for a variety 
of theoretical techniques. The effect of these uncertainties on the total rate measurements 
can be mitigated by measuring partial rates in several regions of phase space. 

In the iriu modes the rates extracted in our chosen q 2 intervals will be largely independent 
of the assumed form factor shapes. In the vector modes, however, the three form factors 
interfere and uncertainties in this interference, particularly for cos 8wi < 0, lead to a form 
factor systematic uncertainty. To reduce this uncertainty, we have included the cos 6 Wi < 
data in the fit, rather than relying on our assumed form factors to extrapolate into this 
region. There is substantial cross feed from the cos 6^/ < region of B — > ptv phase space 
into the reconstructed B — > ntv decays. As a result, the the largest form factor systematic 
uncertainty in the previous CLEO B — > -nlv measurement arose from the B —>■ piv mode. 
Our improved measurements over a larger region of phase space better constrain the B —>■ piv 
rate and reduce this systematic uncertainty by roughly an order of magnitude. 

III. QCD ANOMALY IN B -» rf'^lv 

As mentioned in the introduction, we can check for consistency and probe for new infor- 
mation in quark symmetries by measuring the B — > rjiu and B — > rftv branching fractions. 
Branching fraction measurements for B — > rj'X s [HI H51 HS1 El] yield larger than expected 
results, indicating extra couplings in the B —>■ rf form factor. Measurements at CLEO of 
rf production in T(15*) decays [25] limit the contribution of the rf — > g*g form factor to 
the anomalously high B —>■ r)'X s rate. This unexpected rate may be accounted for by the 
additional QCD Anomaly [26J contribution to the rf axial vector current, which can be 
cleanly measured in semileptonic decays [27] • The QCD Anomaly provides additional gluon 
couplings to singlet states, but no additional couplings to non-singlet states. 

The 577(3) octet of the u, d, and s quarks contain the ir~ , 7r + , K°, K~ , K + mesons as 
well as sums of the tt° with the octet triplet rf = ^(mm + dd — 2ss). The singlet of this 

symmetry is 77° = ^{uu + dd + ss). The physical rj and rf mesons appear to be primarily 

made up of combinations of the 77° singlet state and the rf triplet state. Since the r] mass 
is so much closer to the K masses, it is primarily rf. The rf is mostly the rf singlet, which 
possesses the strong gluon couplings in the axial vector current. A measurement of the ratio 
of B(B — > rflv) to BiB — > r\tv) provides experimental input needed to extract the strength 
of these gluon couplings. 

Once the sizes of these couplings have been determined for a particular model using the 
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semileptonic B — > ??£z/ and -B 
used to check for consistency with the B 



rj'iu decays, the same form factor parameterization can be 
r]'X s decays (see section VII C.) 



IV. EVENT RECONSTRUCTION AND SELECTION 

This study utilizes the full 15.5 fb -1 set of data collected at the T(4S) with the 
CLEO II [28], CLEO II.5 [29] and CLEO III [30] detectors at CESR. The analysis rests 
upon associating the missing energy and momentum in each event with the neutrino four- 
momentum, an approach enabled by the excellent hermeticity and resolution of the CLEO 
detectors. Charged particles are detected over at least 93% of the total solid angle for all 
three detector configurations, and are measured in a solenoidal magnetic field with a mo- 
mentum resolution of 0.6% at 2 GeV/c. Photons and electrons are detected in a CsI(Tl) 
electromagnetic calorimeter that covers 98% of the 4tc solid angle. A typical tt° mass res- 
olution is 6 MeV/c 2 . Unless otherwise noted, all kinematic quantities are measured in the 
laboratory frame. 

Electrons satisfying p > 200 MeV/c are identified over 90% of the solid angle by using the 
ratio of energy deposited in the calorimeter to track momentum in conjunction with specific 
ionization (dE/dx) information from the main drift chamber, shower shape information, and 
the difference between extrapolated track position and the cluster centroid. Depending on 
detector configuration, either time-of-flight (CLEO II/II.5) or Ring Imaging Cerenkov [31] 
(CLEO III) measurements provide additional e ± /K ± separation. 

Particles in the polar angle range | cos 9\ < 0.85 (0.65) that register hits in counters beyond 
five interaction lengths are accepted as signal muons with the CLEO II/II.5 (CLEO III) 
detector configurations. Those with | cos 9\ < 0.71 and hits between three and five interaction 
lengths are used in a multiple-lepton veto, discussed below. 

We restrict signal electron and muon candidates to the momentum interval 1.0 < p < 
2.8 GeV/c and 1.2 < p < 2.8 GeV/c, respectively. The lower limits are determined by 
excessive background contributions for electrons, and by identification efficiency for muons. 
We also require that the lepton candidate tracks have signals from at least 40% of their 
potential drift chamber layers and are consistent with originating from the beam spot. 

Momentum-dependent lepton identification efficiencies and the rates for hadrons being 
misidentified as a lepton ( "fake rates" ) for both electrons and muons are measured in data 
by examining the rate at which the decay products from cleanly identified A — > pn^ , K$ — > 
7T + 7r~, and D* + — > tt + D° — > 7t + K~7t + decays satisfy lepton identification criteria. Within 
the signal electron fiducial and momentum regions, the identification efficiency exceeds 90% 
while the fake rate is about 0.1%. For signal muons above 1.5 GeV/c, the identification 
efficiency also averages above 90%. For muons below 1.5 GeV/c, the five interaction length 
requirement causes the efficiency to fall rapidly to about 30% at the lowest momentum of 
1.2 GeV/c. The muon fake rate is about 1%. The measured efficiencies and fake rates are 
employed in all Monte Carlo (MC) simulations. 

Each charged track that is not identified as an electron or muon is assigned the most 
probable of the pion, kaon or proton mass hypotheses. The probability for each mass 
hypothesis is formed from an identification likelihood, based upon the specific ionization 
measurements in the drift chamber and either the time-of-flight or Cerenkov photon angle 
measurements, and the relative production fractions for pions, kaons and protons at that 
momentum in generic B meson decay. 

To reconstruct the undetected neutrino we associate the neutrino four-momentum p v with 
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the missing four- momentum p m i SS - In the process e + e~ — > T(AS) — > the total energy of 
the beams is imparted to the BB system; at CESR that system is very nearly at rest because 
beam energies are symmetric and the beam crossing angle is small (~ 2 mrad). The missing 
four-momentum in an event is given by p miss = (E miss ,p miss ) = p total - X>char gc d - E Pneutrai, 
where the event four-momentum p to tai is known from the energy and crossing angle of the 
CESR beams and ^charged and p n eutrai are the four-momenta of charged and neutral particles 
explicitly reconstructed by the detector. Charged and neutral particles pass selection criteria 
designed to achieve the best possible |p m i SS | resolution by balancing the efficiency for detecting 
true particles against the rejection of false ones. 

For the charged four-momentum sum X]p c hargcd) optimal selection is achieved with topo- 
logical criteria that minimize multiple-counting resulting from low-momentum tracks that 
curl in the magnetic field, charged particles that decay in flight or interact within the de- 
tector, and spurious tracks. Tracks that are actually segments of a single low transverse 
momentum "curling" particle can be mis-reconstructed as two (or more) separate outgoing 
tracks with opposite charge and roughly equal momentum. These tracks are identified by 
selecting oppositely charged track pairs whose innermost and outermost diametric radii each 
match within 14 cm and whose separation in is within 180° ± 20°. We choose the track 
segment that will best represent the original charged particle based on track quality and 
distance of closest approach information. We employ similar algorithms to identify particles 
that curl more than once, creating three or more track segments. We also identify tracks 
that have scattered or decayed in the drift chamber, causing the original track to end and 
one or more tracks to begin in a new direction. We keep only the track segment with the 
majority of its hits before the interaction point. Spurious tracks are identified by their low 
hit density and/or low number of overall hits, and rejected. 

For the neutral four-momentum sum, ^Pncutrai, clusters resulting from the interactions of 
charged hadrons must be avoided. As a first step, calorimeter showers passing the standard 
CLEO proximity- matching (within 15 cm of a charged track) are eliminated. Optimization 
studies also revealed that all showers under 50 MeV should be eliminated. The processes 
that result in separately reconstructed showers ( "splitoffs" ) within about 25° of a proximity- 
matched shower tend to result in an energy distribution over the 3x3 central array of the 
splitoff shower that "points back" to the core hadronic shower. We combine this pointing 
information with the ratio of energies in the 3 x 3 to 5 x 5 arrays of crystals, whether the 
shower forms a good 7r°, and the MC predictions for relative energy spectra for true photons 
versus splitoff showers to provide an optimal suppression of the splitoff contribution. 

Signal MC events show a |p m i ss | resolution of w 0.1 GeV/c and an -E m i SS resolution of 
~ 0.2 GeV after all analysis cuts. 

To further enhance the association of the missing momentum with an undetected neutrino 
in our final event sample, we require that the missing mass, M^ iss = Ef niss — |p m i SS | 2 , be 
consistent within resolution with a massless neutrino. Specifically, since the M^ iss resolution 
scales as ~ 2_E miss /as miss , for each exclusive B — > X u tv mode we reconstruct, we require 
M^ iss /2E miss , with typical ranges of —0.5 < M? niss /2E miss < 0.3 GeV. The criterion was 
optimized for each decay mode using independent samples of MC simluations for signal and 
background processes. 

Association of the missing four-momentum with the neutrino four-momentum is only 
valid if the event contains no more than one neutrino and all true particles are detected. For 
events with additional missing particles or double-counted particles, the signal modes tend 
not to reconstruct properly, while background processes tend to smear into our sensitive 
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regions. Hence it is worthwhile to reject events in which missing particles are likely. We 
exclude events containing more than one identified lepton, which indicates an increased 
likelihood for multiple neutrinos. 

We require the angle of the missing momentum with respect to the beam axis satisfy 
| cos#| < 0.96. In addition to suppressing events in which particles exited the detector along 
the beam line, this requirement also suppresses events resulting from a collision of radiated 
photons from the beams. These two-photon events are background to BB events and are 
typically characterized by a large missing momentum along the beam axis. 

A nonzero net charge Q indicates at least one missed or doubly counted charged particle. 
The Q = sample offers the greatest purity. For the pseudoscalar modes, use of the \Q\ — 1 
in addition to the Q = sample offers a statistical advantage. Events from the subset in 
which a very soft pion from a D* decay was not reconstructed in particular have only a 
modestly distorted missing momentum. While we retain the \Q\ — 1 sample, we treat it 
separately from Q = sample in reconstruction and fitting so that the statistical power 
of the latter sample, with its better signal to background ratio, is not diluted. For the 
vector modes, which have the poorest signal to background ratios, we require Q = because 
systematic errors associated with the normalization of the larger background in the \Q\ — 1 
sample outweigh the statistical benefits. 

After our selection criteria, the remaining b — > civ background events are dominated by 
events that contain either a Kl meson or an additional neutrino that is roughly collinear 
with the signal neutrino. 

With an estimate of the neutrino four-momentum in hand, we can employ full recon- 
struction of our signal modes. Because the resolution on -E m i SS is so much larger than 
that for |p m iss|, we use (E u ,p u ) = (|p m iss|,Pmiss) for full reconstruction. The neutrino 
combined with the signal charged lepton (£) and hadron (h) should satisfy, within reso- 
lution, the constraints on energy, AE = (E v + Ee + Eh) — -Ebeam ~ 0, and on momentum, 
Mhi v = [£beam ~~ \ a Vu + Pe. + Ph\ 2 }^ ~ Mb, where a, a momentum scaling factor, is chosen 
such that aE u + Ee + Eh — Et, eam = 0. The neutrino momentum resolution dominates the AE 
resolution, so the momentum scaling corrects for the mismeasurement of the magnitude of 
the neutrino momentum in the Mm v calculation. This correction also reduces the correlation 
between Mm v and AE mismeasurement. Uncertainty in the neutrino direction remains as 
the dominant source of smearing in this mass calculation. 

We reconstruct q 2 = M^* = ip v + Pe) 2 for each decay from the reconstructed charged 
lepton four-momentum and the missing momentum. In addition to using the scaled recon- 
structed momentum ap u described above, the direction of the missing momentum is changed 
through the smallest angle consistent with forcing M Mv = M B . This procedure results in a 
q 2 resolution of 0.3 GeV 2 , independent of q 2 . 

A signal charged tt + candidate must have the pion hypothesis as its most probable particle 
ID outcome, and it must not be a daughter in any reconstructed K s decay. A tt° candidate 
must have a 77 mass within 2 standard deviations of the n° mass and an energy greater than 
250 MeV. Each daughter shower must have an energy greater than 30 MeV. Combinatoric 
background levels outstrip efficiency gains below either of these requirements. 

Reconstructed p° — > 7r + 7r~ and p + — > 7r + 7r° candidates are accepted within 285 MeV/c 2 
of the nominal p mass to accommodate the broad p resonance width. The mr mass ranges 
are divided into three 95 MeV/c 2 intervals for the remainder of the reconstruction and 
fitting. We accept uj candidates, reconstructed via its 7r + 7r~7r° decay, within 30 MeV/c 2 
of the nominal uj mass, and, similar to the p modes, the 7r + 7r - 7r° invariant mass region is 



9 



3601206-005 




M(tot) (GeV) 

FIG. 3: The shape of the mass of the p in B — > piv and 7/ — > p^. Limited phase space in rj — > pj 
produces the asymmetry about the nominal p mass in this mode. Histograms are normalized to 
equal area. 



divided into three equal intervals of 10 MeV/c 2 . This division in both modes makes the fit 
sensitive to signal through rough p/u line shape information. Except where noted, plots in 
this paper will include data from only the invariant mass interval centered on the resonance 
mass although data from all three regions are included in the fit. 

We reconstruct 77 in both the 77 and the 7r + 7r~7r° decay modes. For 77, we require the 
reconstructed mass to be within 2 standard deviations (about 26 MeV/c 2 ) of the 77 mass and 
each shower to have cos 9 < 0.81 relative to the beam pipe. Showers within a reconstructed 
7T° are vetoed. For the 77 — > 7r + 7r~7r° submode, the 7r° must be within 2 standard deviations 
of the 7T° mass and statisfy E^o > 225 MeV. To reduce combinatoric backgrounds, the 
invariant mass of the two charged pion tracks must satisfy M nn < M v — M n o + 0.05 GeV/c 2 . 

The 7/ mode is reconstructed in 7/ — > /ry and 7/ — > i]tttt; 77 — > 77. Selection criteria are 
optimized separately based on the event net charge and q 2 where we distinguish between 
q 2 > 10 GeV 2 and q 2 < 10 GeV 2 . Here we specify the requirements for the Q = and 
q 2 < 10 GeV 2 sample, for which we have the greatest sensitivity. For 7/ — > p7, combinatoric 
backgrounds are reduced by requiring 6' 7r7 , the angle between a pion and the photon in the p 
rest frame, to satisfy | cos < 0.925. Photons from n° candidates were vetoed. We accept 
p candidates within the mass range 0.3 to 0.9 GeV/c 2 , where the range is subdivided into 
four roughly equal mass bins. Note that phase space restrictions in 7/ — > p^y distort the p line 
shape at high tttt mass (Figure |3|. The line shape falls rapidly to zero at M n7T = 0.9 GeV/c 2 , 
so the largest mass bin, 0.78 GeV/c 2 — 0.9 GeV/c 2 , is expected to have relatively little 
signal. That mass region is primarily used in the fit to constrain the background level 
at lower mass. The reconstructed 77' mass must be within 2.5 standard deviations of the 
nominal 7/ mass. 

For the 77' — > rj-KTr (7/ — > 77) mode, the 77 daughter photons must satisfy the same criteria 
as the 77 in r]£i>. We also require E v > 300 MeV. Finally, we require the reconstructed 7/ and 
77 masses to satisfy (a 2 , + a 2 ) 1 ^ 2 < 3.75, where cr v (r)') is the number of standard deviations of 
the reconstructed mass away from the nominal mass. 

We minimize multiple candidates per event to simplify statistical interpretation. Within 
each 7r, 77 or 7/ — > i]7i + 7i^ mode, we allow only one candidate within the region 5.175 < 
M Uv < 5.2875 GeV/c 2 and -0.75 < AE < 0.25 GeV used in the fit (Section |v|- Within a 
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FIG. 4: A comparison of |i ca nd • ^roe| for B — > ir£v signal MC (shaded) and off-resonance contin- 
uum data (open). 




FIG. 5: A comparison of R2 for B — » ixlv signal MC (shaded) and off-resonance continuum data 
(open) . 



mode, multiple candidates are resolved by choosing the one with the smallest |Ai£|. For the 
vector and the rf — ► 7r + 7r~7 modes, the \AE\ criterion combined with large combinatorics 
outside the signal region can induce a severe efficiency loss. To mitigate that loss, we select 
the best candidate separately within each of the individual M 7Tn or M 37r ranges. The \AE\ 
criterion induces only very slight peaking in backgrounds, which is modeled by the MC 
generated distributions used to fit the data. 

Backgrounds arise from the e + e~ — > qq and e + e~ — > t + t~ continuum, fake leptons, 
b — > civ, and B — > X u iv modes other than the signal modes. Backgrounds from continuum 
processes make up approximately 75% of the total cross section at the T(49) energy. The 
continuum processes typically produce two collinear jets of hadrons, and the jet-like shape of 
the continuum events can be used to separate those events from the isotropic BB resonance 
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TABLE I: The reconstruction efficiency and fit signal yield for each final state. Both numbers 
are integrated over all phase space, using the nominal form factors for our fit (Section [V]). The 
efficiencies include the branching fractions for the reconstructed modes of the signal hadron. For 
the case of plv and ujIv the efficiency and yield in the central mass bin are reported. Errors on 
the yield are statistical only. 

Decay Mode Efficiency Fit Yield [Events] 



B - 


> T&lv 


4.3% 


179 ± 20 


B - 


■> 7T°£V 


2.6% 


60 ±6 


B - 


r\lv 


1.0% 


14 ±7 


B - 


* p ± lv 


0.8% 


71 ±9 


B - 


-> p°£u 


1.6% 


79 ±10 


B - 


-> u)tv 


0.7% 


34 ±4 


B - 


-> rfiv 


0.5% 


41 ± 12 



events. To do this, we employ a Fisher Discriminant [32J constructed from twelve variables 
that describe event shape. 

The first input variable is |t can d • £roe|, where t can d is the thrust axis obtained using 
our candidate particles (excluding the neutrino), and £roe that for the remaining particles 
in the event. Continuum events typically have hadrons collimated into jets and therefore 
l^cand • £roe| will peak near one (shown in Figure [Zj). For T(45) — ► BB, the B mesons are 
nearly at rest and the thrust directions from the two B decays will be uncorrelated; therefore 
kcand • £roe| will be uniformly distributed. 

We also use the ratio of the second to zeroth Fox- Wolfram moments [33], R2 = H2/ H$. 
The jet structure present in continuum events enhances the second moment and therefore 
the ratio R2 tends to one for jet-like continuum events and to zero for isotropic BB event 
(Figure [5]). 

These two variables provide the main discrimination between BB and continuum events. 
To gain further discriminating power we also consider the angle 9 between the thrust axis of 
the entire event, t ev ent, and the beam axis. For BB events this variable will be randomly dis- 
tributed, while for continuum events the £ ev ent will align with the jet axis and be distributed 
with approximately the 1 + cos 2 9 dependence of the cross section for e + e~ — > qq. 

The remaining nine variables track the momentum flow of the event into nine forward- 
backward cones about the event thrust axis, each spanning 10° in polar angle from the thrust 
axis. We optimize the continuum suppression criterion for each mode (and each subregion of 
phase space for each mode) separately. Due to their jet-like nature, continuum backgrounds 
typically reconstruct with low q 2 and can be largely isolated at q 2 < 2 GeV 2 . The continuum 
suppression has an efficiency of roughly 15% for continuum background and 90% for signal 
B — > iriu events. 

The signal selection efficiencies, averaged over phase space, for the restricted signal region 
range from w 0.5 — 4% and are summarized in Table |lj For cases such as 77 and r/', where we 
only reconstruct a fraction of the decay modes, the branching fractions [37] to the modes 
we reconstruct are included in the efficiency. 
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TABLE II: The q 2 versus cos 8wi intervals for which we extract partial branching fractions for 
each decay mode. 



Mode Index g 2 


rvange [lore v j 


cos 6\yi R-cinge 


iriv 


1 


u < q < z 


1 < COS < 1 


Triv 


o 
Z 


z < q < o 


1 < COS uyyi < 1 


ir£v 


3 


8 < g 2 < 16 


~1 < COS 9\yi < 1 


■kIv 


4 


g 2 > 16 


~1 < COS 6^ < 1 


piv 


1 


< q 2 < 2 


~1 < COS 6^; < 1 


piv 


2 


2 < q 2 < 8 


~1 < COS 6^ < 1 


piv 


3 


8 < q 2 < 16 


< COS 9wi < 1 


piv 


4 


g 2 > 16 


< cos 9wi < 1 


piv 


5 


g 2 > 8 


~ 1 < cos 6 W i < 


Tjtv 




all 


~1 < COS^vi/i < 1 


T)'£l> 




all 


~1 < COS^vi/i < 1 



V. EXTRACTION OF BRANCHING FRACTIONS 
A. Method and Binning 

In order to determine partial branching fractions as a function of q 2 for B — > niv, and 
as a function of both q 2 and cos Owi for £> — > piv we divide the reconstructed candidates 
into the regions defined in Table |TT| The variables AE and Mmv provide sensitivity to signal 
events and therefore in general we bin the data into seven coarse regions spanning the ranges 
5.1750 < M Uv < 5.2875 GeV/c 2 and -0.75 < AE < 0.25 GeV in the AE - M hlv plane. 
We are primarily sensitive to signal within the range defined by 5.2650 GeV/c 2 < M^„ < 
5.2875 GeV/c 2 and -0.15 GeV < AE < 0.25 GeV. The binning choice results from a 
compromise between sensitivity to the signal processes versus reliance on our simulations 
to reproduce in detail the signal and background shapes, which depends critically on the 
modeling of the missing energy and momentum. 

In all of the pseudoscalar modes, we obtain AE versus M Mv distributions separately 
for the \Q\ — 1 and Q = subsamples. This separation allows us to take full statistical 
advantage of the cleaner Q — 0, while obtaining some statistical gain from the \Q\ — 1 
sample. Use of the \Q\ — 1 sample also provides a systematic advantage by reducing our 
sensitivity to simulation of our absolute tracking efficiency. 

We have limited statistics in the 7] and r]' modes; therefore, we extract only total branching 
fractions. However, the backgrounds increase dramatically at large g 2 where the momentum 
of the 7/ ) is low. We choose to separate the reconstructed B — > r/'lv candidates into two bins: 
g 2 greater and less than 10 GeV 2 . Like the division for net charge in the pseudoscalar modes 
this avoids diluting the purer q 2 < 10 GeV 2 sample while still allowing the q 2 > 10 GeV 2 
sample to contribute to the fit. 

Finally, we fit AE versus Mmv distributions for each of the individual M n7T or inter- 
vals discussed above for modes involving a p or u, respectively. The relative yields in these 
regions provide resonance line shape information to help separate signal from backgrounds. 

To extract the branching fraction information, we fit the reconstructed AE versus 
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TABLE III: Summary of the 532 bins used in the nominal fit. The different p mass bins in the 
rj — > /3°7 decay have been categorized as X u decay bins. 



M Mv , AE Q M Xu X u Decay q 2 /6 wl Total 

l^E 7 2 1 1 4 56~~ 

vr°£z/ 7 2 1 1 4 56 

p ± £z/ 7 13 1 5 105 

p°&/ 7 13 1 5 105 

wiv 7 1 3 1 2 42 

r/£z/ 7 2 1 2 1 28 

rj'iv 7 2 1 5 2 140 



distributions for all modes and all subregions of phase space simultaneously using a binned 
maximum likelihood procedure. A summary of the number of bins used in the fit appears 
in Table [TTT1 

The shapes of the AE versus M Uu distributions for signal and background are difficult 
to parameterize analytically. In particular, use of the missing momentum to estimate the 
neutrino momentum produces nontrivial correlations between AE and M^ v even though 
our reconstruction procedure attempts to minimize correlation. We therefore rely on MC 
simulations or data studies to provide the components of the fit (described in detail in the 
following section). To include the finite statistics of these components, our fitter uses the 
method of Barlow and Beeston [31] . 



B. Fit Components and Parameters 

We model the contributions of various processes to the data, outlined below, using a 
combination of MC simulation and independent data samples. All MC samples incorporate 
a full GEANT 3 (35] model of the three generations of the CLEO detector, with event 
samples approximately in the ratio of the number of BB decays from each dataset. The 
simulations also track time-dependent detection efficiencies and resolution. 

Where applicable, our MC simulation incorporates known inclusive and exclusive B decay 
modes as of 2004 [3H]- We correct the generic B model to ensure we can reliably represent 
background processes and signal efficiency losses due to multiple missing particles, primarily 
as a result of Kl particles or charm semileptonic decays within an event. Based on K$ — > 
7r + 7r~ probes of K° production, we find that we must increase the weight of events with 
generated Kl mesons in order to raise the average number of K\ per event by a factor 
of 1.072. We also corrected the inclusive rate and lepton momentum spectrum of charm 
semileptonic decays to agree with the convolution of the inclusive B — > D^*'X spectrum [38] 
with the inclusive charm lepton momentum spectrum from CLEO j3H] • 

Lepton identification efficiency and the rate at which hadrons are misidentified as ( "fake" ) 
leptons are difficult to simulate accurately. In our simulations, we therefore require the 
reconstructed track for a lepton candidate to originate from a true generator level lepton. We 
then apply momentum-dependent electron or muon identification efficiencies obtained from 
independent data studies. We then correct for the efficiency loss that arises in our multiple 
lepton veto due to hadronic fakes. To do so, we randomly veto events using an event-by-event 
probability for an event to contain one or more hadronic fakes. The probability is based on 
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the hadronic content and measured fake rates that are also determined from independent 
data samples. 

We have five categories of components input to our fit. 

1) Signal B — > X u iv Decays: For our nominal fit, we input signal components using 
form factors from the unquenched LQCD calculations by HPQCD [TU] for the iriu modes, 
form factors from the LCSR results of Ball and Zwicky [12] for the p(uS)£v modes, and form 
factors from from the ISGW2 model [21] for the r)iv and rjlv modes. 

We divide the signal MC samples at the generator level into the q 2 and cos 9wi regions 
summarized in Table [TTJ From each of these subsamples we create the 76 AE versus Mm v 
distributions described in Section V A for input to the fit. Thus we input a generalization 
of the full efficiency matrix to the fit, and the crossfeed rates among different regions of 
phase space within a reconstructed mode, as well as among all the different modes, are 
automatically tied to the observed signal yields. The sample sizes ranged from 69 to 300 
times the expected rates. 

Where possible, we use isospin and quark symmetry arguments to combine decay modes 
in order to increase the statistical precision of the extracted decay rates. Within 7r£v, the 
normalization for each of the four q 2 regions for if^l^v floats independently. The neutral 
pion rates are constrained to the charged pion rates to be consistent with isospin symmetry, 
that is T(B° — > 7T~£ + v) = 2T(B + — > 7i £ + u). Similarly the normalization for the five phase 
space regions in p^t v float independently, while the five neutral p rates are constrained 
to be consistent with isospin symmetry. While we reconstruct iotv in only two bins in the 
q 2 1 cos 9wi plane, the signal MC is divided into the same five generated regions as the p 
modes. Based on quark symmetry, we impose the same normalization on the p° and uj 
rates. For our nominal fit we use the 2004 Heavy Flavor Averaging Group [10] values of 
/+-//00 — 1-026 ± 0.034 and t b +/t b o = 1.086 ± 0.017 in applying isospin constraints. 

We fit for the total branching fraction in both the r\lv and the r\tv modes. The branching 
fractions for each 77 and 7/ decay mode are fixed to be consistent with [37]. In rj'iu, the relative 
signal strengths in the two q 2 bins are fixed by the rate predictions of the input form factor. 

In total, the partial rates in each of the phase space bins for B° — > ir~£ + b> and B° — > p~l + v 
in addition to the total B + — > r]t+v and B + — > vj ' i + v are described by eleven free parameters. 

2) b — > c Decays: The dominant b —>■ c background is modeled using a Monte-Carlo 
simulation that incorporates known inclusive and exclusive B branching fractions. To reduce 
systematic sensitivity in simulation of this component, we allow the b — > c normalization 
to float independently for each of the reconstructed modes, net charge bins, and hadronic 
submodes. This procedure introduces fifteen free parameters for a total of 26 free parameters 
in the fit. 

3) Continuum Background: The fit components for the e + e~ — > qq background are 
obtained using data samples collected 60 MeV/c 2 below the T(4S) peak. These samples have 
luminosities of 1/2 (1/3) the total CLEO II/II.5 (CLEO III) on-resonance (BB) luminosity. 
While qq processes dominate this sample, it also accounts for residual two-photon and t + t~ 
backgrounds. In reconstruction, we account for the shift in energy by scaling Mm v by 
-^beam/ -^beam an d) f° r the neutrino rotation in the q 2 calculation, scaling Mb by the inverse 
energy ratio. 

The combination of our continuum suppression and full B reconstruction requirements 
results in low continuum backgrounds overall. As a result, the off-resonance distributions 
have many bins in the Poisson statistics regime. Because the off-resonance data must be 
scaled upward to match the on-resonance luminosity, the Barlow-Beeston extension of the 
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FIG. 6: The Mhi u distribution for B — > 7r^z/ candidates from off resonance data before (points with 
errors) and after (filled histogram) the continuum smoothing procedure. 



likelihood procedure can still result in bias from the very low statistics bins. 

To avoid bias, we adopt a continuum smoothing procedure in which we reanalyze the 
off-resonance data with the explicit continuum suppression criteria removed. This proce- 
dure produces a relatively high-statistics set of M^v — distributions. We correct each 
distribution for potential bias caused by the removal of the continuum suppression. The 
correction is obtained from comparing analysis of an independent mock continuum sam- 
ple with and without the continuum suppression criteria of our analysis. MC provides the 
e + e~ — > qq, t + t~ component with true (generator-level) reconstructed leptons. The fake 
lepton component is obtained from sample of purely hadronic off-resonance data events com- 
bined with measured fake rates. After application of the bias correction, each distribution is 
normalized so that the integrated rate in the Mhtv — A.E region input to the fit is identical 
to the integrated rate obtained from the off-resonance data with the continuum suppression. 
Figure [6] shows one example. 

The off-resonance component is absolutely normalized in the fit based on the ratio of on- 
to off-resonance integrated luminosity with a correction for the energy dependence of the 
e + e~ — > qq cross section. 

4) BB Fake Signal Leptons: Both BB and continuum processes contribute backgrounds 
where a hadron has faked a lepton. Our continuum sample incorporates this contribution, 
but, as noted above, we do not rely on MC for the background component resulting from 
fake leptons. 

We obtain the BB fake component from analysis of on-resonance and off-resonance data 
samples in which no leptons have been identified. For each reconstructed track satisfying our 
electron kinematic and fiducial criteria in each of these events, we analyze the event treating 
that track as the signal electron. If the event with that fake lepton candidate satisfies our 
analysis, then a contribution weighted by the measured electron fake rate for a track of that 
momentum is added. Each track is then similarly treated as a muon candidate. The on- 
and off-resonance samples are absolutely normalized in the fit based on luminosity. The 
off-resonance receives a negative normalization, causing the continuum component to be 
subtracted from the on-resonance data, thereby isolating the BB component. 

5) Other B — ► X u iv Decays: A final background arises from other B — > X u £v decays 
that we are not exclusively reconstructing in the fit. To model this background we use 
a hybrid exclusive-inclusive MC developed and documented by Meyer [12] that combines 
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TABLE IV: A summary of the different components and the corresponding fill style in the subse- 
quent histograms. 





signal 


►♦♦♦♦*♦♦] B — XJv other 




it cross-feed 


| fake lepton 


>//// 

//// 


p/uj cross-feed 


j continuum 




77/77' cross-feed 





ISGW2 predictions of exclusive decays [2T] with the inclusive lepton spectrum predicted 
using HQET by De Fazio and Neubert [43J. The parameters used in the heavy quark 
expansion are constrained by measurements of the B — > X s 7 photon spectrum by CLEO [44J . 
The model generates resonant decays according to ISGW2 plus non-resonant decays such 
that the total generated spectrum matches theoretical predictions. We explicitly remove 
our exclusive decays from this simulation, and the remainder of the sample becomes the 
B — > X° ther £v component of the fit. 

The normalization of this component uses a measurement of the inclusive charmless 
leptonic branching fraction by the BaBar collaboration |35] in the momentum endpoint 
range of 2.2 to 2.6 GeV/c. BaBar finds B(B -> X n e+z/) endpoint = (2.35 ± 0.22) x 10" 4 . 
Roughly 10% of the generated B — ► X° ther £i> spectrum has a lepton in this range. For each 
iteration of the fit, we recalculate the normalization of the "Other B —>■ X u £v" component 
so that the contributions of this component and of the seven signal modes to the endpoint 
region sum to the measured rate. 



C. Verification of Fitting Procedure 

To test the fitter for systematic biases, we performed a bootstrap study using 214 mock 
data samples which were generated by randomly selecting subsets of our MC events. For 
each mock sample, the numbers of events for each signal and background contribution were 
fixed to our efficiency-corrected yields. From fits to the 214 analyzed subsets, we find that 
the fitter reproduces the branching fractions without bias within the statistical uncertainties 
of the procedure (about 7% of the reported statistical error on a given quantity). Further- 
more, the parametric uncertainties reported for the fits agree well with the widths of the 
distributions of the parameters' central values. 

The fit implementation was also cross-checked against an independent implementation 
used in the previous analysis pi] on that analysis' data. 



D. Fit Results 



Table |V| summarizes the branching fractions obtained in the nominal fit for the iriu, p£u, 
r)£v, and rj'ii/ modes. The nominal fit converges with —2 In £ equal to 541 and we note that 
the statistical errors on some of the 532 bins are not Gaussian. Figures |7|- 11 show Mm v 



and AE projections for the nominal fit. These figures are generated by plotting Mmu for 
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FIG. 7: The Mmu and AE projections of the nominal fit, Q = 0, for the summed ir^iv and tPIv 



modes. The fit components are described in Table IV 



candidates with AE in the signal bin and vice versa. Note that for the ir and p modes, the 
binning is finer than that which the fitter uses in order to show detailed peak structure. Raw 
signal yields for the various decay modes, integrated over phase space, are listed in Table [1} 

Figure 12 shows the lepton momentum spectrum for B — > ntv and B — > ptv events in 
the signal AE-M Mu bin. As expected, the signal lepton spectrum for B — > ptv is noticeably 
harder than that for B — > ir£v. The agreement between the data and the sum of the fit 
components gives us confidence in the overall fit quality and our in ability to model the 
lepton momentum distribution of the fit components. 

Figure 13 shows the projections of cos#vw and the invariant mass of the two pions, M n7T , 
used to construct the p candidate for the plv modes. Again, the agreement between data and 
the sum of the fit components is excellent. In the cos 9wi distribution the signal accumulates 
near one and the b —>■ c background is concentrated in the area where cos 8wi is less than 
zero as expected. The M n7T projection shows a clear resonance component and extrapolates 
reliably into regions not used in the fit. The peak in the highest bin is due to B — > D^*Hv 
where the D — ► decay is identified as a p — > tttt decay. 



VI. SYSTEMATIC UNCERTAINTIES 

The systematic uncertainty in this analysis are dominated by effects that couple to ac- 
curate simulation of the missing momentum and energy used to reconstruct the neutrino, 
and to simulation of the selection criteria to isolate events where this association works well. 
Additional uncertainties enter from modeling of the background processes. We explicitly 
evaluate our residual uncertainty due to form factor shapes. 
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modes. The fit components are described in Table IV 



TABLE V: Summary of the phase space subregions and the partial and total branching fraction 
results. The errors shown are statistical and systematic, respectively. 







q 2 [GeV 2 ] 


cos 6 W i 


B [10- 4 ] 






- 2 


-\ - 1 


0.13 ±0.07 ±0.02 






2 - 8 


"1 - 1 


0.27 ±0.08 ±0.03 






8 - 16 


"1 - 1 


0.56 ±0.09 ±0.05 






> 16 


-1 - 1 


0.41 ±0.08 ±0.04 






all phase space 


1.37 ±0.15 ±0.11 


B°- 


* p~t+v 


- 2 


-1 - 1 


0.45 ±0.20 ±0.15 






2 - 8 


-1 - 1 


0.96 ±0.20 ±0.29 






8 - 16 


- 1 


0.75 ±0.16 ±0.14 






> 16 


- 1 


0.35 ±0.07 ±0.05 






> 8 


"I - 


0.42 ±0.18 ±0.31 






all phase space 


2.93 ±0.37 ±0.37 


B°- 


-> 7]£ + V 


all phase space 


0.44 ±0.23 ±0.11 


B°- 


-> rff+u 


all phase space 


2.66 ±0.80 ±0.56 



For most systematic effects, such as simulation of the absolute tracking efficiency, we have 
information from independent studies that limit their size. For most studies, statistics limit 
the sensitivity, so a Gaussian treatment is not unreasonable. To assess the effect of a given 
systematic on this analysis, we completely recreate the fit components with the suspect 
quantity biased to the limit of the study. This typically involves reweighting or biasing of 
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FIG. 9: The M^ v and AE projections of the nominal fit for the summed p^tv and p°£v modes. 



The fit components are described in Table IV 
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FIG. 10: The Mmv an d AE projections of the nominal fit, \Q\ = 0, for the rjiu mode. The fit 



components are described in Table IV 
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FIG. 11: The M^ v and AE projections of the nominal fit, \Q\ = 0, for the rj'iu mode. The fit 



components are described in Table IV The r/ — > p7 and rf — > rjirir submodes have been summed. 



the generated MC samples, which minimizes statistical issues in the determination of the 
shifts. For studies involving use of random numbers (e.g., discarding tracks to decrease the 
tracking efficiency), we repeat the study several times with different seeds to ensure that we 
have an accurate measure. We then refit, and assign the difference from the nominal fit as 
the systematic estimate. Note that we apply the bias to the same samples used to create the 
nominal fit components, which largely eliminates statistical fluctuations in the procedure. 

This procedure will, in principle, cause compensating changes between the efficiencies 
input to the fit and the signal yields before efficiency correction obtained from the fit. De- 
grading resolution, for example, mainly smears signal out of the signal region (decreased 
efficiency), but allows greater latitude for backgrounds to smear into the signal region (in- 
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FIG. 12: Lepton momentum projections of the nominal fit for AQ = 0, n^iv and 7r £i^ (top left), 
p^iv and p ^ (top right), r\tv (bottom left) and rj'iv (bottom right). The data have been summed 



over q 2 , cos9y/h and, where applicable, decay mode. 



creased background, so decreased signal yield). Since the systematic uncertainties can result 
in improper estimation of this cancellation, we will trust any cancellation to only 30% of 
itself, adding that portion back into the assessed uncertainty. 

Because of correlations between yields in different regions of phase space, the system- 
atic uncertainty for the total branching fraction can be smaller than uncertainties for the 
individual partial branching fractions. 



A. Systematic Uncertainties in Neutrino Reconstruction 

The systematic uncertainties associated with neutrino reconstruction efficiency and res- 
olution, summarized in Table VI, dominate the systematic uncertainties in this analysis. 
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cos0 wl M m (GeV/c 2 ) 

FIG. 13: The cos Owl projection of the nominal fit, summed over q 2 , for the fPtv and p^lv 
modes (left). The projection of the two-pion invariant mass of the p candidates summed over q 2 
and cosOyyi (right). Arrows indicate the region included in the fit. The peak at right is due to 
B — ► D(*>£v, with D — * Kir and the K misidentified as a it. 



Track Reconstruction Efficiency: Track reconstruction efficiency agrees in data and sim- 
ulation to within 2.6% (5%) of itself for low momentum tracks (under 250 MeV/c) for 
CLEO II/II.5 (III) and to within 0.5% for high momentum tracks. We randomly discard 
tracks at these levels to assess the systematic effect. We conservatively assume that the 
uncertainties are fully correlated among the three detector configurations. 

Track Momentum Resolution: Based on studies of the reconstructed Kir mass resolution 
in D — > Ktc at the T(4S r ), we increase the deviations between reconstructed and gener- 
ated track momenta by 10% (40%) in CLEO II/II.5 (III). The CLEO II/II.5 value is very 
conservative, while CLEO III value brings data and simulation into proper agreement. 

Shower Reconstruction Efficiency: Shower reconstruction efficiency in data and simula- 
tion agree within 1.6% (1%) of itself for CLEO II/II.5 (III). We therefore randomly discard a 
fraction of the showers consistent with these limits. We again assume full correlation among 
different detector generations when modifying the simulation. 

Shower Resolution: We increase all photon mismeasurements by 10% of themselves. 

Hadronic Showering "Splitoff" Simulation and Rejection: We find that the number of 
splitoff clusters reconstructed from hadronic showers in our simulation deviate from the 
number in data by at most 0.03 per hadron. The energy spectrum agrees well. We therefore 
add additional showers with the observed spectrum at the rate of 0.03 per hadron in each 
event. 

We also systematically distort the neural-net rejection variable, which is based upon 
shower shapes, until the data and simulation clearly disagree. We do this for both showers 
originating from photons and from hadronic shower simulations. 

Particle Identification: Particle identification uncertainties directly affect the missing 
energy resolution. The device-based particle identification efficiencies are understood at 
the level of 5% or better, and our uncertainties as a function of detector and momentum 
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are considered in our systematic evaluation. The overall effect of these uncertainties gets 
reduced by our weighting according to particle production probability 

K\ Production and Energy Deposition: Our Kl production rate correction is modified 
according to the statistical uncertainties of our K s data versus simulation study. A separate 
study of kaons that shower in the Csl calorimeter indicate that the deposited energy is 
simulated to within 20%. 

Secondary Lepton Spectrum: We vary corrections based on the measurements of the inclu- 
sive B — ► D^*'X spectrum and the electron momentum spectrum in inclusive semileptonic 
charm decay according to the envelope of potential b — * c — > stv lepton spectra resulting 
from the uncertainties in the spectral measurements. 



B. Additional Sources of Systematic Error 



The full list of systematic uncertainties are summarized in Table |VTI| and described below. 

Continuum Suppression: We vary the independent parameters in the functions used 
in the continuum smoothing algorithm that parameterize the efficiency of our continuum 
suppression algorithm in the M^ v — AE plane within their uncertainties. 

b — > c: We vary the B — > D£v, B — > D*£u, B — > D**£v, and non-resonant B — > X c tv 
branching fractions according to the uncertainties obtained in Reference [36J. We also apply 
identical variations to the B — > D*lv form factors that are outlined in Reference [4*6] . 

Other B — > X u iv: We vary the heavy quark expansion (HQE) parameters |43j at the heart 
of our hybrid B — > X u iv generator consistent with uncertainties in the CLEO B — > X s ^ 
photon spectrum measurement jH] , which is conservative compared to current world knowl- 
edge. We also probe the hadronization uncertainty by changing the rate for the exclusive 
modes (not including our signal modes) by ±10%, simultaneously adjusting the rate of the 
inclusively generated portion to maintain the same total rate. 

We vary the endpoint rate to which our B — > X u lv component is fixed, within the 
uncertainties of the BaBar endpoint measurement [45] . All uncertainties are combined in 
the summary table, though the HQE and endpoint constraint uncertainties dominate. 

Lepton Identification and Fake Leptons: We have measured our lepton identification 
efficiencies at the 2% level. We also vary the rate at which hadrons fake leptons by ±50%. 
In momentum regions with measurable fake rates, this variation is very conservative. 

7T° Identification: The CLEO III MC simulation overestimates the 7r° efficiency by 4 ±2%. 
Therefore we apply a correction of 0.96 to the CLEO III MC signal decays that contain a 
7T° in the final state. Note that tt finding is only used in the process of reconstructing 
signal decays. To assess the systematic uncertainty, we remove this efficiency correction in 
CLEO III. 

Number of T — > BB Events: Based on luminosity, cross section, and event shape studies 
the error on the number of BB events in CLEO II + II. 5 (III) is taken to be 2% (8%). 
Combining these uncertainties and accounting for the relative luminosities in each, we find 
an uncertainty on the total number of BB events, and therefore the branching fractions, to 
be 3.6%. 

t b+/ t b° an d /+-//oo-' The fixed relative strengths of the charged and neutral rates are 
sensitive to both the lifetime and production ratios of charged and neutral B mesons. Since 
the fit constraints on the charged and neutral ir and p depend on f + _/f 0Q , this systematic 
affects both the total yields from the fit and the conversion factor from the yields to branching 
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fractions. We vary both the lifetime and production fraction ratios independently within 
the uncertainties produced by the Heavy Flavor Averaging Group [%0] . 

Final State Radiation: Final state radiation corrections calculated using PHOTOS [UJ 
are applied to the nominal results. We assign an uncertainty by varying final state radiation 
corrections by 20%. 

Non-Resonant B — > mriv: In the B — ► piv modes there is an additional uncertainty 
from the unknown contribution of non-resonant B — > iririv. Using isospin and angular 
momentum arguments outlined in [1] , we expect the dominant non-resonant irn contribution 
to be in 7r + 7r~ and 7r°7r°. While extracting a true non-resonant 7nriv yield is difficult with 
some knowledge of the form factors and mass distribution, we can limit the amount that a 
p line shape can be projected out of the ir 7r°iv data, and use the isospin relationships to, in 
turn, bound the non-resonant component in 7r + 7r~ iv. Procedurally, we add the reconstructed 
7r°7r°iv mode to the fit, and add one additional parameter that scales the non- resonant 
7t + 7t~£i> and rates. We find we could project out a yield, relative to the charged p 

yield, of —2.2% ± 6.2%, consistent with zero. 



C. Dependence on Form Factors 

The previous CLEO measurement of the B — ► niv and B — > plv rates achieved a very 
minimal dependence on a "priori form factors for B — > niv by measuring independent rates 
in q 2 intervals similar to those in this analysis. Significant dependence on the B — > piv form 
factor, however, remained. Adding the cos 9wi < information in this analysis has greatly 
reduced this residual dependence. 

To evaluate the dependence on the form factors, we choose a set of calculations that span 
the a very conservative range of theoretical results. For each form factor, we re-weight our 
signal MC and refit to assess the impact of the change. We vary the B — > niv and B — > piv 
form factors independently. For each of those, we assign a systematic contribution for the 
rate in each measured phase space region to be 1/2 the difference of the largest and smallest 
rates obtained in that region for the different form factors. The systematic errors due to 
uncertainties in the signal decay form factors are summarized in Table VIII| 



For B — > niv, we consider form factors from the unquenched lattice QCD calculations 
by HPQCD [10] (our nominal form factor), from Ball and Zwicky [UJ, Scora and Isgur 
(ISGW2) |21j . and Feldmann and Kroll [20] . The variation in q 2 dependence of the rate 
from these calculations is illustrated in Figure [TJ 

For B — ► piv, we consider the Ball and Zwicky LCSR calculation (our nominal form 
factors), as well as those of Melikhov and Stech [22], from the quenched lattice calculations 
of UKQCD [23], and from Scora and Isgur (ISGW2) (2JJ. Figure @ shows the the calculated 
rates as a function of q 2 and cos 9wi- Note that the variation between calculations for the 
rate as a function of cos 8 wi is significant, which led to the remaining model dependence of 
the previous analysis. 

Our main sensitivity in the B + — > r]'i + v mode comes from the q 2 < 10 GeV 2 region, 
increasing our sensitivity to this mode's form factor shape. We assess this sensitivity by 
varying the fraction of the B — > rj'iv rate with q 2 < 10 GeV 2 by ±10% of itself. 
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TABLE VIII: The systematic uncertainties associated with the individual form factor uncertainties 
are listed below. The numbers correspond to the q 2 and cos 9y/i intervals that were independently 
reconstructed in the fit. The phase space interval indices are as defined in Table |TT] 
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VII. DISCUSSION OF RESULTS 

A. \V ub \ 

Extraction of \V u b\ from a measured B — > X u £v rate requires theoretical input in the 
form of the form factor(s) governing their decay. The unquenched LQCD calculations for 
the B — > iriv form factor make a minimal number of assumptions. Therefore, our primary 
result for \V u b\ derives from our B — > iriv measurement in combination with LQCD results. 
The p is unstable in unquenched lattice calculations, resulting in a multibody hadronic final 
state that current lattice technology is unable to accommodate. Our primary result for 
B — > piv is a validation of the form factor shapes needed to control the plv backgrounds in 
other B —>■ nlv measurements. 

A determination of \V u b\ from our B — > tt£v results will be most competitive if the 
full range of q 2 can be utilized. However, unquenched LQCD calculations currently cover 
the range q 2 > 16 GeV 2 [10] . The experimental constraints on the form factor shape 
do allow for useful extrapolation of the LQCD results outside this region JUJ HE], but 
techniques for optimal use of both the experimental and LQCD information are still under 
development [19]. Unquenched lattice calculations that utilize a moving B meson frame to 
extend the q 2 range [5U] are also under development, but not yet complete. We therefore 
currently limit ourselves to the range q 2 > 16 GeV 2 , for which the HPQCD collaboration 
finds T(B° -> tt-£+z/; q 2 > 16 GeV 2 ) = \V ub \ 2 (2.07 ± 0.41 ± 0.39 ps" 1 ) [TO]. Using a B° 
meson lifetime of 1.53 x 10~ 12 s [37], we find \V ub \ = (3.6 ± 0.4 ± 0.2+°;°) x 10~ 3 . The errors 
are, in order, statistical, experimental systematic, and theoretical systematic, where the final 
theoretical systematic error is driven by the expected uncertainty in the normalization of 
the form factor as calculated by our chosen LQCD calculation. 



Figure 14 summarizes the current measurements by both Belle [8J and BaBar [61 |7J along 
with the results presented in this analysis for the q 2 > 16 GeV 2 region. There very good 
agreement among the various results and analysis techniques for partial branching fraction 
in this region of phase space which is currently the key experimental input to determinations 
of \V u b\ using unquenched LQCD calculations. 



B. B — > piv Form Factors 



For B — > plv, we perform a single-parameter fit to our measurements over the p phase 
space for the ISGW2 [21], LCSR pi], and quenched LQCD [51] (qLQCD) calculations. The 



28 



3600507-019 



r(B°^7r-| + v;q 2 >16 GeV 2 )/r 



Total 



BaBar [7], 
Untagged, rc* 



This Result, 11°+^ 



Belle [8], |_ 
SL Tagged, ji° 



Belle [8], J 
SL Tagged, jt* 

BaBar [6], ^ 
SL+Hadronic Tagged, xP+n£ 



0.2 0.4 0.6 0.8 
X 10" 4 



FIG. 14: A comparison of the current measurements of the B° — ► partial branching 

fraction in the g 2 > 16 GeV 2 region. The inner (outer) error bars are statistical (statistical and 
systematic). Results from Belle [8] use semleptonic B decay to tag the non-signal B. BaBar has 
published results for an analysis [6] using hadronic and semileptonic B decays to tag the non-signal 
B and also an analysis [7J utilizing loose v reconstruction. Where applicable B° — ► ir~ l + u and 
B + — » ir®l + v results have been averaged assuming isospin symmetry. A weighted average and one 
standard-deviation error band, assuming uncorrelated errors, are shown. 



LCSR results are the results most heavily relied upon for simulation of the pip background 
to The calculation is expected to be valid for the region q < 14 GeV 2 . We therefore 
fit this calculation both to our measurements spanning the full q 2 — cos 8\yi range and to our 
three subregions (see Table |v| that are restricted to q 2 < 16 GeV 2 . The full phase space 
fits test the validity of a given model for use as a background model. The restricted phase 
space fit for LCSR tests the shape without extrapolation outside its region of validity. The 
results are summarized in Table IX All of the probabilities of x 2 exceed 10% for the LCSR 
and LQCD, indicating reasonable agreement between theory and experiment at our current 
precision. The ISGW II model has only a 5% probability. While these constitute the most 
stringent tests to date of the B — ► piv shape predictions, we cannot yet experimentally 
distinguish between the published central LCSR values and the effective shape variations 
used for background systematics in the recent BaBar B — > -niv measurement [7]. 
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TABLE IX: Summary of \ 2 an d number of degrees of freedom (-/Vdof) f° r fits to the B — ► 
differential rate. 



Form Factor 


Phase Space 


x 2 


iVdof 


ISGW II 


full 
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4 


LCSR 


full 
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LCSR 


restricted 


4.3 


2 


LQCD 


full 


4.3 
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C. B{B+ -» ry^+zy) and -» 

We do not see a statistically significant signal in the B — > 77^1/ mode. When B{B + — > 
r]t + v) is forced to zero, —2 In £ increases by 4.8, which corresponds to less than a three 
standard deviation effect. While less statistically significant than reported by the previous 
analysis, the yields in the new relative to the old dataset are consistent within statistical 
errors. 

Forcing B(B + — > ri'C. + v) = 0, — 21n£ increases by 17.9, but we must fold in systematic 
uncertainties to evaluate the significance. Taking the systematic uncertainties as Gaussian, 
we convolute them with our statistical probability distribution function (p.d.f.) inferred 
from the — 21n£ distribution of the fit to model the total p.d.f.. From a toy MC based 
on this p.d.f., we determine the probability for measuring a rate greater than or equal to 
our central value, if the true rate were zero, to be 1.2 x 10~ 3 . Therefore, we interpret the 
significance of our signal to be just over a three standard deviations. 

Using the procedure just described, we find the 90% confidence interval of 1.20 x 10 -4 < 
B{B + — > Vjl^v) < 4.46 x 10~ 4 . For B — > rjiu, we obtain a 90% confidence level upper 
limit of B(B + — > 7]l + v) < 1.01 x 10~ 4 . We also determine the lower limit for the ratio 
B(B + -> r)'l+v)/B(B + -> to be 2.5 at the 90% confidence level. 

Our results allow us to make some model dependent statements regarding the size of the 
extra gluon couplings from the 77° singlet component of the 77 and rj'. Beneke and Neubert use 
the Feldmann, Kroll, and Stech (FKS) mixing scheme [13] in an analysis of B — > 1]'X S [52] 
to model the B — > 77 W form factors as 

F+ =F+ x +F+ w. • (> 

where F B ^ n is the B — ► tx form factor, F smglet is an unknown singlet component, and / (,') 
are constants determined by the FKS mixing scheme. Kim [27] recognized that comparison 
of B decay to r)tv and rjlv could determine the size of the singlet component. 
With the model-dependent parameter 
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5let | 2 0^g 2 







(6) 

where VL^ is the appropriate phase space factor, the B + — > rf''t + u to B + — > n°£ + u branching 
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FIG. 15: The change in the — 21n£ of the fit as a function of the parameter F s . 



ratios may be expressed as 

2B n o 
2B^o 



fa a v , + 27^ a v >b v , y F s + F s b v , 



(7) 



The parameters a (/) and fe w are combinations of FKS mixing parameters, while (3 (,) , 7 w 
and t capture differences in the and _p^ mglet q 2 dependence. Following Kim, we assume 
that the F B ^ and F singlet have the same q dependence, and cover shape differences with 
a systematic uncertainty. 

We impose these branching ratio relationships in our standard fitting procedure and 
allowing F s to float as a free parameter. We find F s = 1.15 ± 0.54 ± 0.38 ± 0.21, where 
the errors listed are, in order, statistical, systematic, and those due to uncertainties related 
to this model, including the relative q 2 dependence. The variation of the likelihood with 
F s is shown in Figure 
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The B -> r^lv branching fractions found with this model are 
consistent with the results of our nominal fit. 

With our F s determination, we can roughly update the B~ — > K~i]' prediction of Benke 
and Neubert by scaling their F smglet (q 2 = 0) = 0.1 (about equivalent to F s ~ 0.137) result, 
for which they found B{B~ -> K'T]') = (56±J2±fJ) x 10" 6 . Experimentally, B(B~ -> 
K~rf) = (71 ±4) x 10" 6 [37]. With our F s , we estimate B(B~ -> K~rf) = (84i|l^) x lO" 6 , 
though a rigorous calculation with F s ~ 1 is needed. This estimate is double the prediction 
with no singlet contribution and agrees well with the experimental value. 



VIII. SUMMARY 

We have measured the branching fractions B (B° — > ir~l + i>) = (1.37±0.15±0.11) x 10 -4 
and B (B° -> p~£ + v) = (2.93 ± 0.37 ± 0.37) x 10" 4 with very little residual dependence of 
either branching fraction on either the B — » 7r or B — > p form factors. Table [V] summarizes 
the branching fraction results for the partial phase space measurements that are summed 
to obtain the measurements integrated over phase space. These results agree well with 
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recent measurements from BaBar [3 El [7] and Belle [8j. The total branching fractions for 
the ixlv and p£v modes are among the most precise current measurements. These results 
indicate that at the level of experimental precision we can probe, the B — > plv form factor 
shapes obtained with LCSR agree with our data, though we cannot test the shape at the 
uncertainty level assumed in the BaBar B — > 7c£u analysis [7j. Using the most recent 
unqenched lattice QCD calculations of the B — > -niv form factor [10J, and measured rate for 
q 2 > 16 GeV 2 region we extract \V u b\ = (3.6 ± 0.4 ± O^i^l) x 10~ 3 , where the uncertainties 
are statistical, experimental, and theoretical, respectively. While this value is competitive 
with other recent determinations of \V u b\, the full strength of this analysis will be more 
effectively realized when a broader q 2 range can be used reliably for the determination of 
\V ub \. 

We find evidence for B + — > i]'£ + u at the three standard deviation level with a branching 
fraction of B(B + -> rf£+u) = (2.66 ± 0.80 ± 0.56) x 10" 4 and a 90% confidence interval 
of 1.20 x 10" 4 < B(B + -> ri'£ + v) < 4.46 x 10~ 4 . The probability that our results for 
B (B + — > r]'£ + v) are consistent with a previous 90% confident limit set by BaBar [17] of 
1.3 x 10~ 4 is approximately 5%. We establish an upper limit for the B + — > 7]£ + v branching 
fraction of B{B + — > 7]£ + v) < 1.01 x 10 -4 (90% confidence level), with systematic uncertain- 
ties included. This is consistent with the BaBar 90% confident upper limit [17] of 1.4 x 10 -4 . 
These results imply a lower limit of B(B + — > 7]'£ + v)/B{B + — > 7]£ + v) > 2.5 (90% confidence 
level). Furthermore, the relative rates indicate a significant form factor contribution from 
the singlet component of the rj', which can bring predictions for BiB — > K^rf) into better 
agreement with experimental measurements. 

These results supersede those previously published by CLEO While some shifts of 
central values are observed, particularly for B° — > p~£v, these results are consistent within 
statistical and systematic errors of those previously presented. We note that retuning of 
the event selection algorithm at low q 2 significantly reduces the statistical and systematic 
correlation between results in this region of phase space. Furthermore, expanding the mea- 
sured B —>■ p£v phase space nearly eliminates the theoretical systematic errors due to the 
B — > p form factor uncertainty that were a significant contribution to the total error on the 
previous result. 

To allow external use of our data in form factor, |T4&|> an d rj and rj' isosinglet studies, we 
provide the full statistical and systematic uncertainty correlation matrices in Appendix A. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. D. Cronin-Hennessy and A. Ryd thank the A. P. Sloan 
Foundation. This work was supported by the National Science Foundation, the U.S. De- 
partment of Energy, and the Natural Sciences and Engineering Research Council of Canada. 



IX. APPENDIX: FULL CORRELATION MATRICES 



Tables [X] and XI provide the correlation matrices for both the statistical and systematic 



errors on the branching fractions outlined in Table |V} Cross-feed among bins produces 
anti-correlations in the results for individual bins; therefore, techniques for extracting \V u b\ 
by averaging over multiple bins in the fit tend to have reduced uncertainty when these 
correlations are taken into account. 
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TABLE X: The statistical correlation matrix for the measured branching fractions. The phase 
space interval indices are as defined in Table [TTJ 



Mode Index 

7T 1 1.000 -0.098 0.001 -0.003 -0.002 0.004 -0.158 0.036 0.008 0.007 0.003 

7i 2 -0.098 1.000 -0.040 0.001 -0.007 0.012 -0.015 -0.195 0.037 -0.003 0.007 

7T 3 0.001 -0.040 1.000 -0.020 -0.025 0.013 0.003 -0.071 -0.173 0.089 -0.105 

7T 4 -0.003 0.001 -0.020 1.000 -0.054 0.018 -0.013 -0.019 -0.015 -0.345 -0.172 

rj -0.002 -0.007 -0.027 -0.054 1.000 -0.008 -0.001 0.002 0.012 -0.058 -0.059 

rf 0.004 0.012 0.013 0.018 -0.008 1.000 -0.010 -0.054 -0.043 0.004 -0.113 

p 1 -0.158 -0.015 0.003 -0.013 -0.001 -0.010 1.000 -0.151 0.031 -0.007 0.048 

p 2 0.036 -0.195 -0.071 -0.019 0.002 -0.054 -0.151 1.000 0.027 0.063 0.102 

p 3 0.008 0.037 -0.173 -0.015 0.012 -0.043 0.031 0.027 1.000 -0.126 0.036 

p 4 0.007 -0.003 0.089 -0.345 -0.058 0.004 -0.007 0.063 -0.126 1.000 -0.267 

p 5 0.003 0.007 -0.105 -0.172 -0.059 -0.113 0.048 0.102 0.036 -0.267 1.000 



TABLE XI: The systematic correlation matrix for the measured branching fractions. The phase 
space interval indices are as defined in Table [TTJ 



Mode Index 

^tT "1 1.000 0.674 0.678 0.383 0.388 -0.033 0.501 0.373 0.337 0.406 -0.257 

7T 2 0.674 1.000 0.754 0.115 0.575 -0.000 0.279 0.012 0.604 0.661 0.246 

7T 3 0.678 0.754 1.000 0.588 0.502 -0.016 0.368 0.258 0.536 0.514 0.056 

7T 4 0.383 0.115 0.588 1.000 0.407 -0.340 0.563 0.704 0.111 -0.099 -0.542 

rj 0.388 0.575 0.502 0.407 1.000 -0.665 0.416 0.406 0.587 0.260 -0.299 

i -0.033 -0.000 -0.016 -0.340 -0.665 1.000 -0.285 -0.477 -0.232 0.228 0.571 

p 1 0.501 0.279 0.368 0.563 0.416 -0.285 1.000 0.852 0.200 0.110 -0.580 

p 2 0.373 0.012 0.258 0.704 0.406 -0.477 0.852 1.000 0.136 -0.159 -0.831 

p 3 0.337 0.604 0.536 0.111 0.587 -0.232 0.200 0.136 1.000 0.737 0.141 

p 4 0.406 0.661 0.514 -0.099 0.260 0.228 0.110 -0.159 0.737 1.000 0.442 

p 5 -0.257 0.246 0.056 -0.542 -0.299 0.571 -0.580 -0.831 0.141 0.442 1.000 
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